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1 Introduction 



Light Front Dynamics is a field theoretical approach which has been successfully applied 
to relativistic composite systems. The two forms of this scheme are: the standard Light 
Front Dynamics (LFD) Q and explicitly Covariant Light Front Dynamics (CLFD) 
While the standard LFD deals with the state vector defined on the plane t + z = 0, this 
plane is defined in CLFD by the invariant equation iv-x = 0, where w is a four- vector 
with = 0. The particular choice of the four-vector to = (1,0,0, —1) turns CLFD into 
standard LFD. 

In this article, we apply CLFD to the calculation, in first order perturbation theory, 
of the QED 7e~ e~ and e~^e~ — > 7 vertices and of the fermion self-energy. We 
will illustrate in details the calculational techniques of CLFD in order to point out its 
differences and similarities with respect to standard LFD and to the Feynman techniques. 

The main difference with calculations in the Feynman approach lies in the fact that 
in CLFD, like in ordinary LFD, all the four-momenta, even in the intermediate states, 
are on their mass shells, whereas the amplitudes may be off-energy shells. Moreover, 
the amplitude may depend, in a well defined manner, on the orientation of the light- 
front plane, i.e., on the four- vector uj. This dependence is an unphysical one for physical 
amplitudes. Thus explicit covariance allows to disentangle clearly the physical amplitudes 
from unphysical ones. In this respect, the CLFD calculations differ from the standard 
LFD ones. 

The study of the perturbative renormalization has already been done in standard 
LFD @-g, resulting in the non-locality of the necessary counterterms. In our covariant 
approach this nonlocality manifests itself only in the terms depending on the orientation of 
the light front plane. These terms can be explicitly removed from the physical amplitude. 
We will show that after their separation, the renormalization of the w-independent part 
of the amplitude is carried out in a very simple way, like in the Feynman approach, 
and does not require any non-local counterterms. Then we find that the on-energy shell 
electromagnetic vertex in CLFD coincides with the on-mass shell Feynman vertex. The 
same is true for the electron self-energy. 

In sect. 2, we start with the calculation of the electron electromagnetic vertex in CLFD: 
the anomalous magnetic moment of the electron, and the renormalized electron charge. In 
sect we apply our formalism to the vertex e"^e~ ^7 for the threshold value of the photon 
momentum = 4m?, and discuss the physical infrared singularity. The renormalized 
electron mass operator is calculated in sect.^. Sect^ contains our concluding remarks. 
Some technical details are given in the appendices |A[and 



2 The electron electromagnetic vertex 

2.1 The anomalous magnetic moment of the electron 

The anomalous magnetic moment of the electron is a simple example of a higher order 
process in QED. Its calculation gives a finite result and does not require renormalization. 
It allows us also to show how to disentangle w-dependent terms in our formalism. 
The spin 1/2 electromagnetic vertex in CLFD has the general form: 

JpiQ) = u{p')Tpu{p) , (1) 
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where q = p' — p. We shall denote = —(^. Due to the explicit covariance of our 
approach, the vertex operator Tp, according to ref.||l^, can be decomposed into: 

Tp = Fi7p + ^cTp^g" ^bA—- \ ) Pp + B^—^p + B^-^Cou^p , (2) 
2m \^^-p [1 + r])m / uj-p [uj-py 

where a^^ = ^(7^7'' — ^'^^f)/2, Co = copj^, r] = Q^/(4m^) and m is the electron mass. The 
electromagnetic vertex (jl]) is gauge invariant since Jpq'' = (with the condition uj-q = 0). 
The possible non-gauge-invariant terms are forbidden by T-invariance. The anomalous 
magnetic moment is the value of F2{Q'^) for = 0. 

The physical form factors Fi and F2 can easily be extracted from the vertex function 
Tp. To this end, we multiply Jp by [vf^' {p')Yu'^ {v)]* ■, [u'^' {p')i(y''" qv / {2m)u'' {p)]* , etc. and 
sum over polarizations. After taking the trace, we obtain the following quantities: 

ci = Tr[0p7^] , C2 = Tr[0piaP''q,]/{2m) , C3 = Tr[Op{u^/io-p - 1/(1 + r])m)]PP , 
C4 = Tr[Op]ujPm/uj-p , C5 = Tr[OpCu]ujPm'^ / (uj-p)'^ , (3) 



Op = {p' + m)Tp{p + m)/(4m2) . (4) 



where 

With the decomposition (^) of Tp, we get a linear system of five equations for Fi,F2, -B1-3 
with the inhomogeneous part determined by 01^5. Solving this system relative to F2, we 
find: 

:[(c3+4c4-2ci)(l+7?)+2(ci+C2)-2(c5 + C4)(l + r/)2] . (5) 



4r/(l -Fr/)- 

In spite of rj in the denominator in eq.(^), there is no singularity at = 0. 

In the usual formulation of LFD on the plane t + z = 0, the form factors of spin 1/2 
systems are found from the plus-component of the current, i.e., in our notation, from the 
contraction of Jp in eq.(||), with ujp. This contraction gets rid of the contributions of -62,3, 
but not the term proportional to Bi. The form factors F[ and F2 inferred in this way are 
thus given by: 

J-uj = u'[FiCj + ^-^ap^ujPq" + 2Bi{Cj - . "^'^ . )]u 
2m [1 + r])m 

= u'iFijp+'^ap.q-'juojf , (6) 

where 

Fi = F, + ^, F^ = F2 + -^B,. (7) 
1 + ?? 1 + r? 

The Bi expression can be found from the above mentioned system of equations, leading 

to: 

^1 = n\ ^ [(C3 + 4C4 - 2ci)(l + + 2(ci + C2) - 4C5(1 + . (8) 
Substituting F2 from (|5|) and Bi from into eq. for F2, we get: 

i^2 = (C5 - C4)/(27?) . (9) 

Of course, in a given order of perturbation theory, both methods for calculating the form 
factors, by eqs.(|5|) and (^), should give the same result. This means we should find Bi = 0. 
We shall see below that it is indeed the case. 
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Figure 1: Electromagnetic vertex of the electron. The dashed line represents the spurion line, 
as explained in details in ref.0 



Let us first calculate the form factor F2, for = 0. The amplitude corresponding to 
the graph of fig. ^ is given by the rules of the graph techniques ||] and has the form: 

/d'T 
u(p)ja(p — k + m)9ioj-{p — k))6({p + ujti — k)"^ — m^) 
Ti — ^0 

X Y{P — k + m)6{uj-(p — k))5((p + L0T2 — k)"^ — m^) ^ — 7uu{p) 

T2 - lO 

X {-g^^'^)e{cO-k)6(e - f.')^ . (10) 

We use the Feynman gauge for the photon propagator. The factor p — k = ki — loti = 
k2 — Cjt2 includes the contact terms —Cjt\ and — a)r2, as explained in ref.Q. For the regu- 
larization of subsequent calculations, we introduced in (|To| ) the photon mass ;U, although 
it is not necessary in the present subsection. 
Integrating over ri, T2 and A;o, we get: 

/ . N 9 /■ u(p)G^u(p) <fk , , 



. (s-m2)2(l-x)2 2efc(27r)3 

where s = {k + ki)'^ = {k + k2)'^, x = uj-k/uj-p and we note: 

GP = --ff,{p -k + m)Y{p -k + m)-fP'. (12) 

The integrands for the scalar functions ci_5 are represented in terms of the scalar 
products between the four-momenta k and uj. The scalar product p-k is given by: 

p-k = + (1 - x)(s - m^)l2 , 

whereas the scalar products w-fc and u;-p always appear in the ratio x, with < x < 1. 

It is convenient to introduce the variable R = k — xp. As usual (see, e.g., ref.p), 
we represent the spatial part of i? as i? = R\\ + where -Ry is parallel to uj and 
is orthogonal to Co. Since, by definition of R, R-co = Rqujq — R^yCu = 0, it follows that 
i?o = and, hence, R^ = —R^ is invariant. In terms of Rj_ and x, the variable s 

writes: 

R'i + ^l^ , R\ + m'^ 



+ -i- — ■ (13) 

\ — X 
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and the integration volume is transformed as: d?k/ek = d? R±dx / x . 
Substituting these expressions into (11), we find (for /i = 0): 

nip)TPnip) = / u(p)G^u(p) f^^'^^ , (14) 

where we denote a = e'^/An. To calculate by eq.(|9p (for r/ ^ 0), we substitute Tp from 
(14) into (0) (for p' = p) and then into expressions (P) for C4 and C5. Calculating the 
traces, we get: 

F^(0) = -^/4A(l-.) (15) 
We thus obtain the well known result for the anomalous magnetic moment of the electron: 

Fm = ^ • (16) 

Now consider the form factor F2 calculated after separation of the w-dependent terms. 
According to (@), it is related to by ^2(0) = F^(0) - 2Si(0). From eq. (|), for = g, 
we find the following expression for Bi: 

- ^ y [iii+m2x2 + ^2(i_^)]2 R^dR^dx , (17) 

which is logarithmically divergent. We regularize it using the Pauli-Villars method, i.e., 
the photon propagator is replaced by: 

111 

(18) 



In the absence of infrared singularity we can put in (^) = 0. Hence, the regularized 
expression for Bi reads: 

B['' = Biiii = {)) - Biiii = K) . 
Integrating it over R±, we get: 



5^' = -— [^dx4- 
Att Jo dx 



m^x^ 



(19) 



After integration over x we get B^^^ = for any value of A. This clearly shows that 
both methods to calculate the anomalous magnetic moment of the electron give the same 
result. 



2.2 The renormalized electron charge 

In order to calculate the radiative correction to the form factor Fi , we have to renormalize 
the charge. The renormalization means that the Lagrangian contains a counter term of 
the form: 

hence, the amplitude is replaced by 

JP Jren ~ ~ -^0 ^ (^0) 
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where 

The renormahzation procedure is described in many textbooks, see for example [11, 12, 
In order to find Zi, one must calculate the amplitude u{p)TPu{p) from the diagram of fig. 
|l[ The value of Zi is in fact just the form factor -Fi(O) determined by this diagram. For 
p = p' the general decomposition (|l|) turns into 



u{p)T'^u{p) = Ziu{p)^^u{p) + Z' u{p)u{p) 



uj-p 



(21) 



where Zi = Fi(0) and Z' = ^2(0) + 53(0). 
From (|2^) the constant Z\ is given by: 



4a; -p 



Tr [bJpT^ip + m)] 



(22) 



The vertex is determined by eq.(^) and is reduced to (O). For regularization purposes, 
we should now keep the photon mass // finite. From eq.(22|) we find: 



a 



(2^)= 



R\ + m?{-2 + 2x + x^ 



R\ + m^x^ + ^^(1 — x) 



-dx . 



(23) 



The subsequent calculation is similar to the calculation given in appendix That is, we 
calculate Zi{fi,L) for a fixed upper limit L of the variable R± in the integral ([2^), take 
the difference Zi(/i, L) — Zi(A, L), take the limit L — > cxo and then calculate the limits 
^ ^ and A — 00. We then obtain: 



Zi (M - 0, A ^ 00) = - + - log ( £^ ) + - log ( 



Svr 2iT 



(24) 



This expression exactly coincides with the expression found in the Feynman formalism 
1 14]. We emphasize that this result for Zi is obtained for the physical part of the full 
vertex (21), after separating out the unphysical term proportional to Z'ujp. The latter 
term can be disregarded, there is no need to calculate it. 



3 Application to the vertex e+e 7 

As a direct application of the preceding calculation, let us now consider the leptonic 
decay width of the positronium. In the Weisskopf-Van Royen limit, the decay width is 
proportional to the elementary vertex e~^e~ — > 7, where the e~^e~ pair originates from the 
positronium wave function with zero relative momentum, i.e. with pg+ = p^- = p. 



3.1 On-energy-shell spin structure 

The on-shell amplitude for the process e~^e~ — > 7 depends on the four-vectors p and uj. 
Its general structure thus reads: 

p^ UJ^TUp 

u{p)M^v{jp) = Au{p)^'^v{jp) + B u{p)Cjv{p) + C-^ yu{p)uv{p) . (25) 
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Here v{p) is the positron spinor. The constant A in ( pq ) is the value of the form factor 
Fi{Q'^) at = Am?. One can also construct the structure aP^^ujp/uj-p, but it is not 
independent, since: 



im 



-u(p)a^^ujgv(p) = u(p)^^v(p) — -^u{p)u)v{p) . 

OJ-p LO-p 

Multiplying (^) on the left by u{p) and on the right by v{p) and summarizing over 
polarizations, we get the factors J2x'^ip)^ip) = (p + w-)) J2x'^ip)^ip) = ip ~ We 
introduce therefore the quantity: 



MP = {p + m)MP{p - m) , 
and calculate the following traces: 



Ti = 

12 = 

Ts = 



1 



1 
1 



Tr 
Tr 

Tr 



MPQj 



P^ 
uj-p 



{^■Pf 



(-3A + C)/2, 
{B + C)/2 , 

= {A + B)/2 . 



So we can find out the coefficients which determine the amplitude (|25|): 

A = T2-T1-T3, 
B = Ti - Ts + 3T3 , 
C = -Ti + 3T2 - 3T3 . 



(26) 



(27) 



(28) 



3.2 The physical ampUtude 

The off-energy-shell amplitude depicted in fig.§ is given by the rules of the graph tech- 
niques |0] and writes: 



u{p)M^v{p') 



u{p)'jf^{k + m)0{u!-k)5{k'^ — m? 



'(2^ 



(29) 



\2 2 

^ m 



X Y(^m-{Q- k)) 0{LoiQ - k))5 [{Q - k + wra) 
X -f.v{p'){-9''n^ ((P - ^r' + con - kf - e{ujip - k)) 



dT2 



T2 - iO 



Ti — iO 

where Q = q — ujt. Note that the fermion and antifermion propagators in LFD differ 



from each other. The propagator (k + m) in (29) corresponds to the electron, whereas the 
propagator (^m — {Q — k)^ corresponds to the positron. 

The factor m — {Q — k) = m — {ki — CJT2) incorporates the difference ki — C0T2 and, 
therefore, takes into account the contact term — a)r2. We consider in this section the case 
q2 _ 4j^2 j-gigvant for the decay width of the positronium, so that p' = p, t' = and so 
the two graphes corresponding to the two different time orderings should give the same 
contribution MP = Mf + = 2Mf . 
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Figure 2: Light-front time ordered graph for e~^e ^ 7. A similar diagram with the opposite 
time-ordering for the photon exchange should be added. 



After integration over ri and T2, the amplitude (29) is given by: 

u(j))Mfv{p) = J u{p)^^{k + m)'^^ {m — {Q — k)^ ^^v{p) 

e{ujip - k))e{ujiQ - k))e{uj-k) d^k 



where 
and 



LU-Q 

si2 -Q^ = 2{uj-Q)t2, S123 -Q^ = 2(wQ)ti , 



(30) 



X]^ 1 Xf^ 

S123 = {k + k2+p) = \ h — 

Xj^ Xj^-^ Xp' 

= S^ + ^k±if! + 2m^ (31) 
Xk 1/2 - Xk 

Like in the previous section, we define above the variables Ri = l — xiQ with xi = uj-I/uj-Q, 
where I is either k, ki or p'. At the threshold = 4m^, we have Rp'± = and Xp/ = 1/2 
in the variable si23- We thus find: 

MP = -2-^ / dPR^ / -f ^— -— , (32) 

{2^Y J ^Jo (si2-4m2)(l-x) (si23-4m2)(l/2-x)2x ' ^ ^ 

where: 

QP = {p + m)-f^{k + m)Y (m - (Q - fc)) 7''(p - m) . (33) 
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The amphtude is connected to MP by eg. (pq) . The factor 2 in ( |3^ ) results from 
the sum of the two amplitudes Mi and M2. In order to find the coefficients A,B,C 
which determine the amplitude (P5|), we substitute M^ into eqs. (H), (HI), regularize the 
expression by the Pauli-Villars prescription: 

A ^ Aif,) - A{A) , 

(and similarly for B, C), and take the limits /i — > 0, A — > 00. The details of the calculation 
are given in the appendix The final expression for A is then: 

am 7a a ^ ( u?\ a ^ ( h?\ 
A = — + — log +— log — . 34 



This expression exactly coincides with that calculated in the Feynman formalism [14|. 
The integral for B{^) converges and does not depend on /x: 

B{^.) = . (35) 

The amplitude regularized a la Pauli-Villars is determined by the difference Breg = B{ii) — 
B(A). It is therefore zero. 



From (34) and ( p4D we recover the renormalized amplitude 



„ am 2a , , 

Aren = A-Z= , 36 

H vr 



which coincides with the result found in the Feynman approach |14]. It contains the term 



corresponding to an infrared singularity. In the next section we will show that this 
singularity is physical but does not contribute to the relativistic correction to the decay 
width of the positronium. 

The calculation of C gives a divergent result even after a single Pauli-Villars regular- 
ization. However, C is the coefficient in front of the term uj^ulju which is proportional 
to ujP. Like Z' in eq. (^ij), it is separated out in eq. ( p5| ) and does not contribute to the 
observable amplitude, determined by A. 



3.3 Infrared singularity 

The infrared singularity in (^) results from the Coulomb interaction between the electron 
and the positron. It manifests itself in the low momentum limit in the loop of fig. 2. 

To calculate this limit, it is sufficient to take the nonrelativistic limit for all the four 
momenta. It means that we make in the numerator of eq.(|30|) the following replacements: 

u{p)^^{k -|- m) — > u{p){l + 7o)w. and {m — {Q — k)j 7^w(p) — m(l — 7o)w(p). 

We used the fact that in this equation only the matrix 7^ with /u = contributes. After 
this substitution, the integral ( pO| ) converges and its calculation, in the /Li — > limit, gives 

atn 

2u(j>)Mfv{p) = u{p)Yv{p)- 

This Coulomb contribution has to be removed from the total radiative corrections (|3^ ) 
since it is already taken into account in the calculation of the positronium wave function. 
After this, we reproduce the well known radiative correction 

This result enables us to calculate the relativistic radiative corrections beyond the 
Weisskopf-Van Royen approximation. This is done for instance in ref.[^] for the calcula- 
tion of the leptonic decay width of the J/ip- 
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4 The electron self-energy 



Since our formulation of LFD is explicitly covariant, we are able to follow very closely the 
standard procedure of renormalization of the fermion self-energy in perturbation theory. 
The self energy diagram is shown in Fig. 3. As already done for the electromagnetic vertex, 
we can write down immediately the general spin structure of the self energy. It is very 
simple and given by: 

= ^i(p2) + i?i(p2)Z + Ciip^)^ . (37) 
m 

The coefficients Ai, Bi, Ci are scalar functions of only. Here p = pi — wti is the total 
momentum entering the diagram, pi is the external fermion momentum, with pi = m?, 
and LOTi is the external spurion momentum. The off-energy-shell mass operator S(p) 
depends on the value p^ = — 2{uj-p)ti. 

Like in the above calculations, we subtract from S(p) the w-dependent structure, 
introducing: 

tip) = S(p) - Ci(p2)^ = A,ip^) + B,ip^)^ . (38) 

m 

The standard procedure of renormalization of Feynman diagrams relies on two countert- 
erms: the mass counterterm Jm^ and the wave function renormalization proportional to 



Z2 |11|. They can be incorporated explicitly in the Hamiltonian. 



Following 1 13], the renormalized self energy S^(p) is defined as the part of S(p) which 
is of second order in the variable (p — m). Without loss of generality, we can rewrite S(p) 
in the form: 

tip) = Ao + ip- m)Bo + Sfi(p) . (39) 

Here Aq,Bq are constants (they do not depend on p^) , and S_r(p) is the renormalized 
self-energy written as: 

^dp) = iP - m)^Mip) , (40) 
where the matrix A4ip) can be represented as: 

Mip) = a + ip + m)b . (41) 




Figure 3: The electron self energy graph 
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The renormalized fermion propagator now reads: 



p — m — {p — mYM.{p) 



It has the same pole and the same residue at ^ = m as the free propagator of the physical 
fermion. 

The explicit calculation of the renormalized fermion self energy is now straightforward. 
According to the rules of CLFD, the electron self energy shown on fig. |3| has the form: 



x9{lo ■ {p - k))5{{p + UJT - kf 



LVT + m)Y{-9fMu) 
a, dr di^k 



iO (27r)3 



{2-Kf 



Am-2k + 2lot d'^R±dx 



■ p^ 



2x1 



(42) 



As indicated in the previous section, we introduce the photon mass fi for infrared regu- 
larization. The term lot i n (p2|) contributes to Ci(p'^)a) only and can be omitted in the 
calculation of T,{p). In eq.(^2|), r = (s — p'^)/{2u)-p) and 



X 



1 



hp 



R^: 1 



— XT) 

2 



with X = uj-k/u> ■ p, R = k — xp, and the phase-space volume is given by d^k/sk = 
d'^R±dx/x. 

Now, starting from eq.(p^) we can calculate and the scalar coefficients a and 

b in (41). Knowing S(p) from eq.(H2]), we calculate the coefficients Ai,Bi in ( |37| ) and 
find S(p) by eq.(|38|) . Comparing ( pq ) with (|39|), we express Ao,Bo through Ai,Bi for 
j)2 = Using again the representation ( |39D for S(p), we finally obtain the functions 
a, b, which determine the self-energy (|40|), through Ai{p'^), Bi[p'^) and -^o- The details 
of the calculation are given in appendix For the functions a and b we find: 



a 



1 



47rm (1 — p) 

a 



2Tim?p 



2(1 -p) 



2-3p 
1-P 

2 



p^+^p-A \ 
P H ; log p +1 + log 

l-p / 



(43) 



where 



m 



P 



With these expressions for a and b the renormalized mass operator writes: 
^r{p) = (P — m)^AI(p) = {p — m)^[a + (p -|- m)6] . 



(44) 



It exactly coincides with the standard result given for instance in refs.|pT|, [l^. Note that 
in higher order calculations, the w-dependent term Ci{p^)Co may give an uj independent 
contribution, when Ti{p) enters as a part of a more complex diagram. It can also be 
renormalized. The calculation of C]^^ is given in appendix B. 
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4.1 The antifermion self-energy 

Since the forms of the fermion and antifermion propagators are different in LFD (they 
contain (j5 + m) for fermion and —{p — m) for antifermion), the form of the self-energy 
is also different. However, there is no need to repeat the calculation: the antifermion 
self-energy can be found from the fermion one. 
We represent similarly to (|37|): 

= A^[p^) + B^{p')^ + Ci{p')(j 
m 

= 'Ao + Bo{p + m) + Ci {p')Cj + , (45) 

where Aq, Bq are constants. The renormalized antifermion self-energy S/j(p) is represented 
in the form similar to (^^: 

'^r{p) = (j5 + m)'^[a — {p — m)b]. (46) 

One can easily see that 

a = a b = b , (47) 

where a,b are given by eq. (|4^ . Indeed, from the comparison of the fermion and an- 
tifermion propagators we see that the antifermion self energy S/j(p) can be obtained from 
the fermion one T,ji{p) by: 

T,ii{p, m) = —Yjr{p, —m) = —{p + m)^[a(— m) + {p — m)b{—m)\ . (48) 

Comparing ( ^6| ) with ( ^8| ) and taking into account that a{—m) = —a{m), b{—m) = b{m) 
(see the explicit expressions (|4^)) we reproduce eqs. (|4^) and (|47|) . 

5 Conclusion 

The understanding of perturbative renormalization in QED is an unavoidable step before 
studying more subtle systems like QCD. While this perturbative renormalization is now 
a text-book section for the standard formulation of field theory using Feynman graph 
techniques, it is not as well understood in Light-Front Quantization. The main reason 
being the difficulty to exhibit the covariant structure of the electromagnetic vertex and 
electron self-energy since standard LFD breaks explicitly covariance. 

We have shown in this study that the covariant formulation of LFD is a powerful tool 
in order to make the link between LFD and Feynman approaches. The explicit covariance 
of our formulation enables us to exhibit the relativistic structure of the electromagnetic 
vertex in QED, as well as the electron self-energy. We are thus in a position to extract, 
after renormalization, the finite physical contribution from the infinite amplitude. To do 
that, we have to know the dependence of the operators on the orientation, of the light 
front. This is trivial in CLFD. In the standard formulation of LFD, this dependence can 
not always be disentangled from the physical part of the amplitude. This w-dependent 
contribution is responsible for the non-locality of the counter terms needed to renormalize 
the infinite amplitude in LFD g-g. 

The finite physical amplitude we found in our approach for the electromagnetic vertex, 
the electron self-energy, and the e^e~ ^7 amplitude agree thus with the standard text- 
book results. 
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We emphasize that in order to reproduce these results, a covariant regularization of 
divergences (Pauli-Villars in the present study) is important. The attempt to regularize 
the integrals by a cut-off in the variable i?j_,for instance, allows of course to work with 
finite integrals. It gives finite, but however wrong, renormalized results. 

The QED vertex e~^e~ — > 7, corrected by a color factor, coincides with the QCD vertex 
q~^q~ 7. This result is applied in ||l5| for the calculation of the relativistic radiative 
correction to the J/ip leptonic decay width. 

The question of non-perturbative renormalization, for scalar particles, will be adressed 



in a forthcoming publication [16| 



Acknowledgements 



One of the author (V.A.K.) is sincerely grateful for the warm hospitality of the Laboratoire 
de Physique Corpusculaire, Universite Blaise Pascal in Clermont-Ferrand, where this work 
was performed. This work was partially supported by the grant No. 99-02-17263 of the 
Russian Fund for Basic Researches. 



A Calculation of A, B and C 



As explained in sect. p.2| , we substitute MP into eqs.(|27|), (|28D in order to find the coeffi- 
cients A,B,C determining the amplitude (p5|). We thus find: 



{A,B,C) 
with: 



(27r)3 



1/2 



(a, b, c) 



dx 



(si2 - 4m2)(l - x) (si23 - 4m2)(l/2 - x) 2x 



t2 — tl — t'S 
1 



^— \ Tr [QPCb] ^-Tr [QP-fp] - Tr 



16m^ I " ' uo-p 
R]_{1 - 2x)+m^{l + Ax'^) 



MP 



tl-t2 + 3t3 



1 



Rj_{l - 4x) + m2(l - 2xf{l + 2x) 



-tl + 3t2 - 3t3 



1 



R\ + 2m^Ri]_{l - ^x^) + m^(l - 4x^)2 



4m^x^ 

To calculate the traces ([50|), we need the following scalar products: 

A;-Q = 2m2 (1 -x)(si2 -4m2), k-p = k-Q/2, p-Q = Q"^ /2 = 2rr? . 
Substituting @ into (||), we find: 

r-1/2 



Bit,) 



Svra 
(2^ 

Svra 
"(2^ 



d'^R± 



d^Ri 





R'j_{l - 2x) +m'^{l + Ax'^) 




Rl + m^{l - 2x)2 




R'{ + m^{l -2x)'^ + 2fi^x 



1/2 





Rl{l 


- 4x) + 


m^(\ - 2x)2(l 2x) 




Rl+m'^{l 


- 2xf 




R\ + m^{\-2xf ^2y?x 



(49) 



(50) 



(51) 



■dx , (52) 
dx ,(53) 
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(2^ 



1/2 





'R\ + 2m'^R\{l- 


8x2) + m''(l 


- 4x2)2 




2m?x 


R\ + m?{l - 2x)2 




R\ + m?{l 


- 2x)2 + 2ij?x 



■dx . 
(54) 



The integral (54) for C, which is the coefficient in front of the structure proportional to 
ujp, diverges quadratically at R_\_ —>■ oo and logarithmically at x = 0. The integral (|K 
for A logarithmically diverges at R± 
asymptotic expression: 



— > oo. The integral (53) for B at R± oo has the 

B{fi) oc / — ^ / (1 - 4x)dx . 
Jo R± Jo 

Since the integral over x is zero, B(fj,) is finite. 

The regularization of A{iJ,) proceeds as follows. The integral over x in A{fj,) can be 
done analytically. In the Pauli-Villars regularization scheme, we should take the difference 
A{fj,) — A{A) and calculate the convergent integral over R±. Equivalently, but technically 
easier, we calculate A(/x,L) with the cutoff L in the variable R±, take the difference 
A{fi, L) — ^(A, L) and then take the limit L — > oo. The result is analytic, but lengthy. In 
the limits ^ ^ and A — > oo we obtain eq.(34), which coincides with the result calculated 
in the Feynman formalism [14]. 



B Calculation of a, 6 and C[^^ 

According to eqs.(^0[) and (|4l|), the self-energy S_r(p) is determined by the scalar functions 
a and b. As we will see below, a and b are determined by the coefficients Ai{p'^) and i?i(p2) 

2 

m . 



(55) 

(56) 
(57) 



in the general decomposition (|3^) and by their combinations in the limit p 
From (|37| ) we find the coefficients Ai(p'^) and Bi{p'^): 

1 TY) 

= -rr[S(p)], Br = —TT[T.{p)(b] . 
4 uj-p 

Substituting here eq.(]4^) for S(p) we get: 
Ai{p') = ^ 
Blip") -- 



irdlf', dx 



7r2 J R\ + (1 — x)m2 + 2;[/i2 + (1 — x)p^] ' 
am I -ndR^, xdx 



27r2 J R\ + {1- x)m2 + x[/j,2 + (1 - x)p^\ 
These integrals diverge logarithmically. 



Comparing (|37D with (|39D and taking into account (|46|), we find: 

P 



A,{p') + Bi{p' 



m 



Ao + {p- m)Bo + (p- mYM{p) . 



From here we can express the constants Aq and Bo through Ai and By. 



Ao 
Bo 



1 

4m 



■Tr 



A,ip^) + B^{p^)^]{p + m) 
m . 



4m (p2 — m"^) 



Tr 



Aiip^) + Bi{p 



2x P 



m 



= Ai(m2) + Bi{m'^ 
Ao]ip + m)2 



(58) 

(59) 
(60) 
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We thus obtain: 

^0 
Bq 



am 



(2 - x)iTdR^, dx 



27r2 J Rl + {l- x)2m2 + xfx^ ' 



Q 



x[Rj_ - (3 - 4x + x^)m^ + x^i^]-KdR\dx 



[R\ + (1 - a;)2m2 + x/i2]2 

These integrals also diverge logarithmically. 

From (|5^), and taking into account eq.(41) for fA{p), we get: 

{p - mf{a + {p + m)b) = Ai{p^) + Bi{p^)— - Aq - {p - m)Bo 

m 

This allows finding out a and b: 
1 



(61) 
(62) 



-Tr 



Ai{p^) + Bi{p^) - Ao 



Ai(p2) + Si(p2)Z _Ao-{p- m)Bo ) (p + m)p 
m 



p2 _ 



4p2(p2 _ 77^2)2 



Tr 



.2n P 



Ai{p') + Bi{p^)^ - Ao - {p - m)Bo {p + mYp 
m 



2m(^i(p2) - ^0) ^ (p2 + rn2)5^(p2) 



p2 — 77^2 



(p2 — r7T,2)2 m(p2 — r?T,2)2 

Substituting here the above expressions for Ai,Bi,Ao and Bq, we get: 

x{2 -3x + x'^)dR\dx 



am 



a 



27r J [R\ + m'^{l-xf][R\ + m^{l-x){l-{l- p)x) 
x2(l - x)[R\ - m2(3 - 4x + x2)](ii22^dx 



(63) 



(64) 



(65) 



2tT J [R\ + m2(l - x)2 + //2a;]2[^2^ _^ 7,^2(1 _ 2:)(1 - (1 - p)x) + fi'^x] 

We omitted in a the photon mass /u, since that integral has no infrared divergence, and 
introduced the notation: p = (m^ —p'^)/m'^. Integrating over i?^ and x and keeping in b 
the leading term in \og{iJ? /m^) only, we obtain eqs.(^3[). 

One can similarly calculate the coefficient Ci determining the w-dependent part of 
S(p). It is given by: 



Ci(/ 



Auj-p 



-Tr 



T.{p) p 



P^Cj 
uj-p 



[2R\ + m2(2 - 3(1 - p)x'^)]'^dR\dx 
Air'^LO-p J + m2(l — x)(l — (1 — p)x) + fi'^x]x 



(66) 



It is quadratically divergent in the variable R± and is logarithmically divergent at = 0, 
despite the finite photon mass /U. Note that the standard Pauli-Villars regularization is 
not enough to make it finite. 

We renormalize this scalar function in the standard way |11]: 

Clip' 



(p2) = Clip') - Ci(m^) - (/ - m') 



dp^ 
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After that it becomes finite: 

am'^p'^ f [Rl{2 - 5x) +m'^{2-5x + 3x^) - 3i^'^x^]x{l - x)T:dR\dx 



Qren 



2n ^ P f 

' 47r2(a;-p) J 
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